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Abstract 

We analyse the multiplicity distributions of charged particles at 
Tevatron (pp) and LHC (pp) energies in the framework of the inde- 
pendent pair parton interactions (IPPI) model. It is shown that the 
number of soft pair parton interactions (and therefore the density of 
the partonic medium) is large and increases with energy. The mean 
multiplicity at each parton interaction grows also with energy. This 
growth depends on the width of the rapidity window. The similar 
conclusions are obtained in the multiladder exchange model (QGSM). 

1 Introduction 

Following Feynman, two colliding strongly-interacting high-energy objects 
(protons, nuclei etc) are considered as sets of pointlike partons. Parton in- 
teractions determine the outcome of collisions with many particles produced. 
They can be either soft, with small transferred momenta, or hard, at large pt- 
Most secondary particles are created in soft processes. Usually their char- 
acteristics are described within some phenomenological models with their 
Monte Carlo (MC) implementations. Particles (or jets) with high transverse 
momentum are more rare and their production is described by perturbative 
QCD albeit at additional assumptions. 

One of the general and widely debated (for reviews see [H El [3] ) charac- 
teristics reflecting mostly the properties of soft processes is the multiplicity 
distribution. At energies of tens of GeV the experimental data are well fitted 
by the negative binomial distribution (NBD) first introduced in jl]. However 
at higher energies this fit by a single NBD becomes inadequate. A shoul- 
der structure appears at high multiplicities. It is tempting to ascribe it to 
multiple interactions of pairs of colliding partons. This was done in several 
models at the expense of introducing new parameters. We have proposed 
[5] the most economic model with minimum adjustable parameters which we 
call the independent pair parton interactions (IPPI) model. In fact, we show 
that there are only two such parameters. Those are the maximum number 
of active parton pairs at a given collision energy and the average multiplicity 
of particles produced in the collision of a single pair. At comparatively low 
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energies one pair is active and it leads to NBD of produced particles which 
fits experimental data. To deal with the shoulders at higher energies it is 
assumed that the number of active parton pairs from colliding particles in- 
creases with energy. This is in accordance with increase of gluon densities 
at low shares of momentum x. Interaction of each additional pair results 
in NBD for its products with the same parameters as for a single pair since 
the interactions are independent. One is actually tempted to recognize that 
this assumption could be valid at asymptotically high energies where energy 
conservation becomes unimportant and the number of colliding partons is 
infinitely large. We'll try to use it at energies of Tevatron and LHC. 

It is the well known property of NBD easily seen from its generating 
function (see, e.g., [2]) that the convolutions of NBDs again lead to this 
distribution (see [5J). NBD is characterized by two parameters - mean mul- 
tiplicity m and dispersion k. With a maximum number of active pairs (jmax) 
at a given energy added, we are left with these three values. However, we 
show that one of these parameters, namely the dispersion, can be excluded 
due to some special property of the IPPI model. Thus, only two of them 
should be used for fits of particle multiplicity distributions at any energy. 
To compare, there are seven parameters (including j ma x) in the multiladder 
exchange or quark-gluon string model (QGSM) [6] used in [7j which are to 
be fitted by additional energy dependent sets of experimental data on total 
and elastic scattering cross sections. Even the simple fit with two NBDs asks 
for six parameters to be used if there are no additional constraints. 

Therefore, the IPPI model is the most economic one concerning the ad- 
justable parameters and can serve as a first approximation for estimates of 
the global characteristics of multiple parton interactions. The IPPI model 
does not imply that there are no correlations between particles. They are 
intrinsic in each binary collision because of assumed NBD and in their convo- 
lution. Surely, further correlations between these interacting pairs of partons, 
of both dynamical and kinematical origin, can be introduced as it is done, for 
example, in the multiladder exchange model [6j which asks for new parame- 
ters. Also, there are more rare processes where some pairs interact strongly 
and are scattered at large transverse momenta (see,e.g., [SI IS])- That would 
lead to 2-, 4- and more-jets (or 66,7 etc) production. However, being inter- 
esting by itself, the latter aspect is out of reach of the present study. 

Let us stress that the knowledge of the number of active parton pairs 
provides information about the density of the hadronic matter formed in 
collisions and about its evolution with energy. We shall see that this density 
is high and grows with energy increase. These findings favor the enlarged 
role of collective effects at LHC energies. 
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2 Application of the IPPI model 

The main equation of the IPPI model obtained in [5] is 

Jmax 

P(n;m,k) = ^ WjP NB D(n;jm,jk). (1) 
i=i 

It states that the probability of the n-particle production channel is defined 
by the sum of NBDs with shifted maxima (jm) and larger widths (jk) for 
independent parton collisions weighted by their probabilities Wj. The total 
number of particles n is equal to the sum of particles produced in all active 
pairs collisions. At asymptotically high energies the probability for j pairs 
of independent interactions Wj is the product of j probabilities for one pair 
so that the normalization condition 

Jmax Jmax 

E *i= E w i = 1 ( 2 ) 

j=l 3=1 

determines w\ if j max is known at a given energy. In fact, the value of wi 
ranges between 1 at low energies (for j max = 1) and 0.5 at asymptotics where 
jmax tends to infinity. Thus, all values Wj are calculated from Eq. (j2J) if jmax 
is defined. We show them in Table I for 6-10 active pairs which happen to 
be important at TeV energies. 

The IPPI model predicts new special features of moments of the mul- 
tiplicity distribution which impose some constraints on the parameters and 
allow to get rid of one of them, namely k. The factorial moments of the 
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distribution (jTJ are 

F q = ^P(n)n(n - !)...(„ -q+l)= ( j)' = /,(*) (f )' 

(3) 

with 

Jmax T^i jh -U Jmax 

m = E ^ ^ = * E + ( 4 ) 

Herefrom one gets the relation 

(5) 

This relation states that the righthand side with the definite ratio of in- 
dependent functions should be independent of q for some value of k. It 
opens the way to the combined fit of experimental multiplicity distributions 
with the requirement of g-independence of m (see Appendix 1). 

This requirement is very strong and can be satisfied at some special values 
of Wj unknown to us. It would be too naive to expect it to be precisely 
satisfied in a simplified model. Nevertheless, we try to find such fits of 
experimental multiplicity distributions by the IPPI model which minimize 
the decline of m as a function of q from constancy. After doing this, we 
get the important information about the maximum number of active parton 
pairs at a given energy and its evolution with energy. It provides the clear 
insight into the dynamics of soft interactions and properties of the hadronic 
medium formed during the collision. 

Beside the factorial moments F q we use the cumulants 

9-1 fa- II' 

K < = F « ~ E -jr 1 } -^K q _ r F r (6) 

and their ratios 

H q = K q /F q (7) 

which possess some specific oscillating behaviour (see [2j [TUl E] ) . This is the 
complementary (and sometimes more sensitive!) approach. 

We have used both direct fits of multiplicity distributions and the fits of 
H q moments at energies of Tevatron 1.8 TeV (for pp interactions) [12] and 
LHC 0.9, 2.36 and 7 TeV (for pp interactions) [33]. The detailed description of 
our procedure is presented in Appendix 1. The extrapolation of the obtained 
results admits predictions at 14 TeV which we also show below. 
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Figure 1: 



The fit of the multiplicity distribution at 1.8 TeV was first done in our 
paper [5]. Now using the improved procedure described in Appendix 1 we 
confirmed the stability of previous fits and their parameters as seen from good 
fit in Fig. [TJ The values of fit parameters are k = 4.42 and m = 12.944 ±0.04 
which practically coincide with those in [5]. The corresponding fit of H q is 
presented in Fig. E] (k — 4.36, m = 12.90 ± 0.05). The values of k and m 
are close in both approaches. The independence of m on q is satisfied within 
the limits less than 1 % as seen in Fig. El Within these limits the value of 
k does not play a role and can be considered fixed by constancy of m at a 
given energy. 

The maximum number of the pair parton interactions is j ma x(l-8 TeV)~4. 

The same procedure was used for the LHC pp-data at 0.9, 2.36 and 7 
TeV. The experimental distributions are available for limited pseudorapidity 
intervals. Namely they have been used by us in distinction to the above 
fit at Tevatron energy where the data were given for the total pseudora- 
pidity window with some extrapolation to the fragmentation region. Beside 
some arbitrariness of such extrapolation, the multiplicity distribution for the 
full phase space is influenced by energy-momentum and charge conservation. 
These constraints are less important in restricted domains. Therefore, the 
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Figure 2: 
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Figure 4: 



distributions can be expected to be more sensitive to the underlying dynam- 
ics. 

The parameter m shows now what part of the multiplicity produced in 
a single pair parton interaction reaches the analysed pseudorapidity inter- 
val. Since intervals \rj\ <2.4 cover the region of almost flat pseudorapidity 
distribution one would expect that the ratios of corresponding values of m 
are approximately equal to the ratios of the intervals themselves. On the 
contrary, in the fragmentation region \q\ >2.4 the distribution drops down. 
Thus the values of m for 2.4 must be close to the parameters obtained from 
extrapolated Tevatron data. All these features as well as the explicit en- 
ergy dependence of m are clearly seen in Fig. H] where some results for lower 
energies from our paper [5] are also shown. 

In the same Figure we plot the numbers of active parton pairs from our 
fits for each pseudorapidity interval near corresponding points. They are 
quite stable for a fixed energy and rise from 4 at 0.9 TeV to 6 at 7 TeV. 
The corresponding contributions to the total mean multiplicity (equal to 
31.4 at 7 TeV for the interval \q\ <2.4) of these 6 pairs interactions are 
8.3+8.3+6.3+4.2+2.7+1.6. 

The achieved good fits of the shapes of multiplicity distributions with 
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Figure 5: 



above parameters assure that the energy dependence of mean multiplicity 
and higher moments are also well reproduced. 

We have compared our conclusions with those which one would obtain 
from the multiladder exchange model (QGSM) [6] used in [7J. The lower 
mean multiplicity of a single pair parton interaction and, correspondingly, 
somewhat larger numbers of ladders are allowed there due to the wider spread 
of probabilities Wj. For example, at 7 TeV they would be 13.6 with 8 active 
pairs of partons. For these fits we have chosen the values of 7 adjustable 
parameters exactly equal to those used in [7J and can not say how sensitive to 
their variations are the results. The overall fit of the multiplicity distribution 
for \rj\ <2.4 according to the multiladder model is somewhat better than in 
the IPPI model (see Fig. [5]). Corresponding % 2 /dof are 62/127 and 131/127. 

The qualitative conclusions about the energy increase of the mean mul- 
tiplicity at a single parton interaction (see Fig. [6]) and the number of such 
interactions (see Fig. [JJ are strongly supported in both approaches. 

In Fig. [8] we plot the predictions at 14 TeV for \rj\ <2.4 in both models. 
They are obtained by using the extrapolated values of m shown in Fig. [6j 
The predicted numbers of active parton pairs at 14 TeV range from 7 in IPPI 
to 10 in QGSM. In general, we see that the difference between these models 
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Figure 8: 

becomes noticeable only at energies as high as 7 TeV even though it is still 
not very well pronounced in Fig. [HJ 

We conclude that the high density partonic medium is formed not only 
in heavy-ion collisions but also at high energy pp-interactions. The num- 
ber of active parton pairs taking part in soft pp interactions increases with 
energy reaching values 6 - 8 at 7 TeV and could be 7 - 10 at 14 TeV, i.e 
their density is higher at higher energies. Therefore the theoretical account 
of multiple parton interactions is absolutely necessary at high energies. The 
mean number of particles created at a single interaction of the pair of par- 
tons also increases with energy. This is quite natural because the structure 
functions are modified correspondingly. The collective effects should become 
more pronounced at LHC energies. 

3 Appendix 1 

The experimental data on multiplicity distributions were obtained at Teva- 
tron (for pp) and LHC (for pp) in various pseudorapidity intervals. Their mo- 
ments can be computed and compared with predictions of the IPPI model. 
To get smaller numerical values at the intermediate stages of computing, we 
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actually deal with factorial and cumulant moments normalized by (n) q (q— 1)!. 



F q n i-, n/r — 1 

^ = <n>*(g-l)! = ^ H M (8) 

= -Fg - £ —fcq-r- (9) 

r=l r 

The values of i7 ? are not changed by this normalization. 

To estimate the errors in moments induced by the experimental error bars 
of P n one usually assumes that the latter ones are independent for different 
n and gets 



: A ^) 2 = E (ff ) ( Ap n) 2 = £(n(n - 1) . . . (n - g + 1)) 2 (AP„) 2 , 



(10) 



(A^) 2 = E(^^^J (ap„) 2 , (ii) 

where dK q /dP n is obtained by differentiating fl§]). 

Moreover, one should include the requirement of g-independence of m 
in the fitting procedure minimizing the decline of values given by Eq. (j5J) 
from a constant when calculations are done within some model. Then one 
minimizes the following functions 



/plPPIfh. ™\ _ p \ 2 Imax 

mm / P (fc, m) = mm }_l — + ^ 

fc,m k,m n \ AP n J 2 



'm(q, k) — 
v Am(q,k) t 



q=2 

(12) 

. /flT ^ . ^/ gf^-^ y ^ m(g,fc)-(m(fc)) V 

mm fS fc = mm > — - — — + ; — — . 13 

k JH K J k ^ 2 \ AH q ) \ Am(q,k) J K } 

Here m(q, k) is computed according to (jHJ) in the IPPI model. As m in 
Eq. f|12p is one of the parameters to be minimized it does not contain experi- 
mental errors (in contrast to (m(k)) in Eq. ffTB"]) ) and the denominator in the 
second sum in Eq. f lT2"j) can be computed as 

. , . m(q, k)AF 
Am(^) = ^ q -. (14) 

The case of Eq. fll3p asks for a more complicated formula which is not 
shown here. Let us note that all the sums are not normalized to the number 
of particles and the number of the moments since it enlarges the relative 
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weight of the constancy of m while our goal is to get the best fit of the 
multiplicity distributions at the satisfactory fit of additional conditions. 

Using these errors one can find out the parameters of the IPPI model in 
the twofold way: either by the direct fit of multiplicity distributions P n , or 
by fits of computed values of H q . We show the results of both approaches. 
The latter one is more transparent and the error bars are easily visualized. 

There is another problem which becomes important. Namely, one must 
decide what is the maximum rank (the value of q m ax) to be used in Eq. 
(j!2p . To answer this question we introduce a simple criterium related to the 
experimental error of AH q . We compute the dispersion of AH q according to 
the standard expression 



Df = (Atf 2 ), - (AH)l (15) 



where 



and 



(AH) q = (AH) q -! + ^(Atf, - (AH) q .{) 

with (AH)! = {AH 2 } 1 = by definition. 

Afterwards we don't consider the moments starting from such rank q that 
satisfies the condition 



AH q - (AH) g 



> 1. (16) 



And in any case we don't consider ranks larger than 16. The maximum 
number of the pair parton collisions j max is then chosen automatically so that 
to minimize the value in Eq. (1121) . In that respect our procedure improves 
the approach of [5], where this parameter was chosen just among two nearest 
possibilities. After that one can apply this procedure in attempts to fit 
experimental data. At 7 TeV we introduce the cut at q max =ll (see Figj9]). 
Namely these results are described in the main content of the paper. 



4 Appendix 2 

The g-independence of m is a crucial test of the quality of fits. Let us just 
briefly mention that we tried to apply the minimization procedure without 
any condition of constancy of m, i.e. with no last terms in Eqs. (TP2"]) and 
([13]) . It showed slightly different (albeit within the limits of less than 10 
per cents) values of parameters k and m if the Tevatron data at 1.8 TeV 
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Figure 9: 

were used. Visually, these fits looked quite satisfactory (see Fig. [10]) both for 
multiplicity distributions and for H q . However the conditon m(q, /c)=const 
is not well fulfilled. Minimization of probabilities lead to values /c=3.79 and 
m=12.59, while minimization of H q gives fc=4.08 and m=12.828. Fig. [TT1 
demonstrates the g-dependence of m(q,k) (j^J) for A;=3.79 and k=4.08 (the 
straight lines show the average values). It is much more noticeble than that 
in Fig. [3j Therefore the last fit discussed in the paper is preferable. 
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